Universal quantum gates based on a pair of orthogonal cyclic states: Application to 

NMR systems 
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We propose an experimentally feasible scheme to achieve quantum computation based on a pair 
of orthogonal cyclic states. In this scheme, quantum gates can be implemented based on the total 
phase accumulated in cyclic evolutions. In particular, geometric quantum computation may be 
achieved by eliminating the dynamic phase accumulated in the whole evolution. Therefore, both 
dynamic and geometric operations for quantum computation are workable in the present theory. 
Physical implementation of this set of gates is designed for NMR systems. Also interestingly, we 
show that a set of universal geometric quantum gates in NMR systems may be realized in one cycle 
by simply choosing specific parameters of the external rotating magnetic fields. In addition, we 
demonstrate explicitly a multi-loop method to remove the dynamic phase in geometric quantum 
gates. 
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I. INTRODUCTION 

Building a practical quantum computer with a large 
number of qubits has recently attracted much attention. 
For realization of a universal quantum computer, there 
are certain minimum requirements: the storage of quan- 
tum information in a set of two-level systems-qubits, the 
processing of this information using quantum gates, and 
a mean of final readout [1]. So far, a number of systems 
has been proposed as potentially viable quantum com- 
puter models, including trapped ions [2], cavity quantum 
electrodynamics [3], nuclear magnetic rcsoncc(NMR) [4], 
and low-capacitance Josephson juctions [5-8], etc.. 

An essential requirement in quantum computation is 
to maintain quantum coherence in a computing system, 
since the coherent interference pattern between the mul- 
titude of superpositions is necessary for taking advantage 
of quantum parallelism. However, the coupling of a quan- 
tum system to its environment leads to the so-called deco- 
herence process, in which encoded quantum information 
is lost to the environment. The error rates of the indi- 
vidual gate operations should be less than 10~ 4 to assure 
that the quantum computer works fault-tolerantly [9] . To 
accomplish the required precision, the decoherence time 
of the system has to be much longer than the operation 
time required for computing. How to suppress the infa- 
mous decoherence effects is one main task for quantum 
computing. 

One of schemes to correct the errors caused by decoher- 
ence is quantum error-correcting codes [10-12], through 
which originally encoded information can be recovered 
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by suitable encondings and measurements of qubits. An 
alternative approach to avoid decoherence has been pro- 
posed in Refs. [13,14], where decoherence-free states have 
been used as qubits. The decoherence-free space is a 
subspace which is inherently immune to unwanted noise. 
In addition, evolution of the system must not drive the 
state out of the decoherence-free space. So far, all of 
these strategies require extra-physical resources, such as 
additional manipulations or encoding one logical qubit 
by several ancillary physical qubits. 

Another attractive strategy for fault-tolerant quantum 
computation is based on a topological idea [15], where 
gate operations depend only on global features of the con- 
trol process, and are therefore largely insensitive to local 
inaccuracies and fluctuations. A significant advance in 
this direction is made by the so-called geometric quan- 
tum computation [16]. In this kind of scheme, a universal 
set of quantum gates may be realized by pure geomet- 
ric phases, which depend only on the geometry of the 
path executed [17-19], and therefore provides a possibil- 
ity to perform quantum gate operations by an intrinsi- 
cally fault-tolerant way [16,20,21]. 

Several basic ideas of adiabatic geometric quan- 
tum computation by using NMR [20], superconduct- 
ing nanocircuits [8], trapped ions [21], or semiconductor 
nanostructures [22] were proposed. However, the adia- 
batic evolution appears to be quite special, and thus the 
nonadiabatic correction on the phase shift may need to 
be considered in some realistic systems as it may play 
a significant role in a whole process [18,23]. A seri- 
ous disadvantange of the adiabatic approach is that the 
evolution time must be much longer than the typical op- 
eration time To of the qubit system, while the evolution 
must be completed within the decoherence time, which 
leads to an intrinsical time limitation on the operation 
of quantum gate. Therefore, a generalization to nona- 
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diabatic cases is valuable and important in controlling 
quantum gates. 

Recently, nonadiabatic geometric quantum computa- 
tion has been proposed theoretically [24,25], and de- 
tection of the conditional nonadiabatic geometric phase 
shifts for quantum gates using NMR is also experimen- 
tally reported [26]. Nevertheless a systematic study of 
this topic, especially the appication to NMR systems, is 
still highly desirable. In this paper, we propose an exper- 
imentally feasible nonadiabatic scheme to achieve a uni- 
versal set of quantum gates [27] based on a pair of orthog- 
onal cyclic states. In this scheme, quantum gates may 
be implemented based on either the total phase accumu- 
lated in the cyclic evolution or the geometric Aharonov- 
Anandan (AA) phase shift [18] after eliminating the dy- 
namic phase. Therefore, quantum computation imple- 
mented by dynamic and geometric operations can be 
unified in the present theory. In addition, physical im- 
plementation of this set of gates is designed in detail for 
NMR systems, in which the qubits considered are spin- 
1/2 particles in the presence of a magnetic field rotating 
uniformly around a fixed axis. Although the Schrodingcr 
equation of this system was solved exactly long time ago 
[28], and the AA phase was previously obtained explic- 
itly in Refs. [23,29], we here derive all phase shifts explic- 
itly and show that they may be applicable in achieving a 
universal set of logical gates. Moreover, the nonadiabatic 
geometric computation may be experimentally achieved 
just by simply choosing specific controllable parameters, 
with the cyclic states as a pair of dark states [24] . 

The paper is organized as follows. In Sec. II, we discuss 
general aspects of the geometric phase and cyclic evolu- 
tion, and then present a theory applicable for achieving a 
universal set of quantum gates based on a pair of orthog- 
onal cyclic states. In Sec. Ill, the theory is applied to a 
viable NMR quantum computer. The paper ends with a 
brief summary. 



II. IMPLEMENTATION OF QUANTUM GATES 
WITH A PAIR OF ORTHOGONAL CYCLIC 
STATES 

For universal quantum computation, we need to 
achieve two kinds of noncommutablc singlc-qubit gates 
and one nontrivial two-qubit gate [27]. Thus we here 
consider only two-qubit systems. A general Hamiltonian 
for two qubits may be expressed as 



i=i 



(1) 



where (k = 1, 2) arc the Pauli operators for qubit k, 
B[ k \t) are local (real or fictitious) magnetic fields act- 



ing on fc-th qubit, and Ji represents the strength of the 
interaction between two qubits. 



A. Cyclic evolution and geometric phases 

Before the design of quantum gates, we present first 
general aspects of the cyclic evolution in a qubit sys- 
tem. A cyclic evolution is referred to as that the state 
of the system returns to its original state after evolution. 
Mathematically, a normalized state \4>{t)) is cyclic in the 
interval [0, r] if and only if 
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IV>(0)>, 



with 7 being a real number. The total phase 7 acquired 
in the evolution would contain both geometric and dy- 
namic components, denoted as "f g and 7^, respectively. 
Usually, not all states take cyclic evolutions. A sufficient 
but not necessary condition for cyclic evolution is that 
the initial state is a nondegenerate energy eigenstate of 
a cyclic Hamiltonian which changes adiabatically. 

We first focus on the cyclic evolution in a qubit system. 
At this stage, only one term 



H = -fjB- a /2 



(2) 



needs to be considered. Here B denotes the total mag- 
netic field felt by the qubit, which may include a real 
external field as well as an effective magnetic field in- 
duced by the interaction between different qubits. The 
Hamiltonian H is chosen to go through a cyclic evolution 
with period r in the parameter space {B}. 

We here present how to calculate all kinds of phase 
shifts. The dynamic phase may be directly calculated 
from its definition given by 



Id 



H\m)dt, 



(3) 



while the geometric part is not easy to derive. We first 
address a method to calculate adiabatic Berry's phase. 
By adiabatically changing H around a circuit in {B}, 
the eigenstate will accumulate an adiabatic Berry's phase 
7b = =fO s /2, where the signs =F depend on whether the 
system is in the eigenstate aligned with or against the 
field, and Q s is the solid angle subtended by the magnetic 
field at the degeneracy B = [17]. U s can be derived as 



f 

Jo 



B x dtB y — B y dtB x 

mBT+W) ' 



(4) 



under the condition of a closed trajectory with B(r) = 
B(0) [25]. 

Since the requirement of the adiabatic evolution could 
be stringent, a generalization to nonadiabatic case is 
more desirable. The generalization of adiabatic Berry's 
phase to a nonadiabatic cyclic evolution was introduced 
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in Ref. [18], where a general geometric phase 7 9 = 7 — 7^ 
is defined as 

i 9 =ij\m\^- t m)), (5) 

Rere_\i>(t)) = e- if ^\tp(t)) with /(r) - /(0) = 7, leading 
to IVK 7 ")) = IV'(O)}- The A A phase can be regarded as 
a geometric phase associated with a closed curve in the 
projective Hilbert space, and approaches Berry's phase in 
the adiabatic limit. The A A formulation applies regard- 
less of the Hamiltonian H being cyclic or adiabatic; it 
depends only on the cyclic evolution of the system itself. 

Normally, Eq. (5) is not directly used to calcu- 
late the geometric phase accumulated in a cyclic evo- 
lution. We here present an alternative method to cal- 
culate the non-adiabatic geometric phases. This ap- 
proach is more convenient for qubit systems discussed 
in this paper. For a spin-1/2 particle in the pres- 
ence of an arbitrary magnetic field, the nonadiabatic 
cyclic A A phase is just the solid angle determined by 
the evolution curve in the projective Hilbert space-a 
unit sphere S 2 . Any two-component 'spin' state — 
{e-^l 2 cos(0/2), e l¥,/2 sin(0/2)] T may be mapped into a 
unit vector n = (sin^cos^s, sin sin 95, cos#) in the pro- 
jective Hilbert space via the relation n = (ip\ a\ ip), where 
T represents the transposition of matrix. By changing 
the magnetic field, the AA phase is given by 



1 — cos 9) dip, 



(6) 



where C is along the actual evolution curve on S 2 , and 
is determined by the equation 



d t n(t) = -fiB(t) x n(t)/h. 



(7) 



This 7 9 phase recovers Berry's phase in adiabatic evolu- 
tion [23]. The cyclic evolution implies that n(i) under- 
goes a closed path in the projective Hilbert space. 

We consider a process, in which a pair of orthogonal 
states \tp±) can evolve cyclically starting from \ip±(0)). 
A pair of orthogonal state may be parametrized as 



and 



IV>-> 



-iip/2 



sinf e*v/ 2 



-sinfe-*/ 2 
cos § e'v/ 2 



(8) 



(9) 



Denoting n±(i) = (ip±(t) \ <r\ ip±(t)), it is straightforward 
to find that n + (i) = — n_(t) by using Eqs.(8) and (9). 
For a cyclic evolution, \tp±(r)) — e n± \ip±(0)). Besides, 
we have an important relation: 7+ = —7-. This is be- 
cause the dynamic phase 



ld+ = -lJ^ E+(t)dt = ~J -E_{t)dt 



= ~7d- 



with 



E±(t) = (il>±{t)\H\i>±) = -/xn±(t) • B(t), 



and the geometric phase 7 S (— n(0)) = — 7 s (n(0)) at any 
time if the two initial states correspond to ±n(0) [23]. 
By taking into account the cyclic condition for \ip±), we 
finally have 



U(T)\ip±) = exp(±i7)|^±), 



(10) 



where U(t) is the evolution operator. Hereafter we de- 
note 7,7 3 , and -fd as the phases for \tp + ) for brevity. 



B. Quantum computation 

We now show how to realize a universal set of quantum 
gates based on either the total phases or the geometric 
AA phases accumulated in cyclic evolutions. 



1. Quantum logical gates 

A quantum logical gate is a unitary operator U act- 
ing on the states of a certain set of qubits, that is, U 
may be referred to as a quantum gate if |Vw) = U\4>in) 
with \ipi n ) being the input state and \ip ut) being the 
output state. The space of all the possible input and 
output states makes up the Hilbert space of states for 
the quantum computer. If Ti is the Hilbert space of 
a single qubit, and is a given basis state for the 
ith qubit, then a basis vector \ip) for the states of the 
quantum register is a tensor product of qubit states 
W) = IV'O ® 1^2} ® • • • ® IVn) G W®"- U N is an N- 
qubit gate when \ipi n ) & Ti.® N . Nevertheless, we need 
not to implement all Ui (1 < I < N) but only a universal 
set of gates U u . A set of gates U u is called universal if 
any unitary action £/; can be decomposed into a product 
of successive gates in U u . It is shown that two noncom- 
mutable one-qubit (single-qubit) gates and one nontriv- 
ial two-qubit gate consist of a universal set of gates [27]. 
This universality is very useful in practice since it allows 
us to focus only on how to construct a universal set of 
gates. 

We first construct the single-qubit gates by assuming 
that a pair of orthogonal states \tjj±) can evolve cyclically. 
We write an arbitrary input state as \ipi n ) = a +\ip+) + 
a-\tp-) with a± = (ip±\ip), and express the two cyclic 
initial states as \tp + ) = cosx/2|0)+sinx/2|l) and = 
— sinx/2|0) +cosx/2|l), where |0) and |1) constitute the 
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computational basis for the qubit. Using Eq.(lO), the 
output state at time r is found to be [25] 

Wout) = U S «{ X ,imn), 

where 

TTH „, . / e* 7 cos 2 $ + e^* 7 sin 2 ^ isinYsin7 

U q (x,l) = • • • 2 i-y • 2 X , - i-y 2 

y i sin x sin 7 e 7 sin f + e 7 cos 

(11) 

For this gate, there exists a relation 

[u sq (x,i)V = [u sq (x,-i)}, 

where the adjoint operation t corresponds to transposi- 
tion and complex conjugation of matrix. Thus the impor- 
tant f operator for a unitary operation U may be exper- 
imentally achieved by the operation U with the inverse 
sign of the phase factor. 

It is straightforward to verify that two operations 
U^'ixiili) an d ^ (X2,72) are noncommutable as long 
as 

sin 7i sin7 2 sin(x2 - Xi) ^ 0. (12) 

Since two kinds of noncommutable operations constitute 
a universal set of single-qubit gates, we achieve the uni- 
versal single-qubit gates by choosing xi X2 + jir for 
any nontrivial phases 71 and 72 (71,2 7^ jn), where j is 
an integer. For example, the phase- flip gate U^^x) — 
exp(— 2i7i|l)(l|) (up to an irrelevant overall phase) is ac- 
complished at x = 0; the gate U^i^) = cxp(z72<Xr) is 
obtained at x = 7r/2, which produces a spin flip (NOT- 
operation) when 72 = ir/2 and an equal- weight super- 
position of spin states when 72 = 7r/4. U^ q 2 are two 
well-known single-qubit gates. 

In terms of the computational basis {|00),|01),|10), 
1 11}}, where the first (second) number represent the state 
in the controlled (target) qubit, the unitary operator to 
describe the two-qubit gate is given by [8,25] 

U tq — diag(U ( 7 °, x °), J7( 7 i iX i)), (13) 

under the condition that the control qubit is far away 
from the resonance condition for the operation of the 
target qubit. Here 7^ (x S ) represents the total phase 
(the cyclic initial state ) of the target qubit when the 
control qubit is in state S(= 0, 1). Following Ref. [27], we 
find that unitary operator (13) is a nontrivial two-qubit 
gate if and only if 7 1 ^ 7 or x 1 ^ x° {mod 2ir). For 
example, 

U tq 0rfl) = dia 5 (e i7 V* 7 V 7 V V ), (14) 

when x 1 — X° — 0; this gate was proposed to be achieved 
in the adiabatic case in the charge qubit [8]. Combining 



gate (14) with single-qubit operations we are able to per- 
form a gate described by 

U C N = [I® U s %*/A, tt/2)] U tq ^ /2) [I ® U s %k/A, 71/2)]* 
= diag(I,ia x ), (15) 

with / as a 2 x 2 unit matrix. This gate is equiv- 
alent to the controllcd-NOT (CNOT, which is defined 
as \m)\n) — > |m)|m ® n), where denotes the addi- 
tion modulo 2.) gate up to an overall phase factor 
for the target qubit. On the other hand, t/(o 7 i) b c ~ 
come the controlled-phase (CPHASE, which is defined as 
|m)|n) — ► \m) exp(imn<fi)\ri) .) gate by removing a overall 
phase for the target qubit. 

An alternative practical method to achieve the con- 
trolled two-qubit gate is also available under certain con- 
ditions. Denoting the Hamiltonian of the target qubit 
as H t , we may produce H t — by choosing certain pa- 
rameters of H t when the controlled qubit is in the state 
|0), while H t is able to realize a required gate when the 
controlled qubit is in state |1). Then the gate in this case 
is given by 

U* = diag(I,U(%x)), (16) 

where 7 is the total phase accumulated in the evolution 
when the controlled qubit is in state |1). Gate (16) cor- 
responds to gate (13) for 7 = and xo = tt/2. Ucn in 
Eq.(15) may be directly derived when 7 = x = 7r /2- 

So far, we have demonstrated that all elements of quan- 
tum computation may be achievable by using a pair of 
orthogonal cyclic states. 

2. Geometric quantum gates 

The quantum gates U described in Eqs. (11), (13) [or 
(16)] may be divided into two categories: one is referred 
to as a geometric gate if the phase in U is a pure geo- 
metric evolution operator (7^ = 0) [30], and the other 
is referred to as a dynamic gate as long as there exists a 
nonzero phase induced from dynamic origin (i.e., 7^ ^ 0). 
Geometric quantum computation demands that logical 
gates in computing are realized by using geometric phase 
shifts, so that it may have the built-in fault-tolerant ad- 
vantage due to the fact that the geometric phases depend 
only on some global geometric features. 

A key point in geometric quantum computation is to 
remove the dynamic phase. We here address two meth- 
ods [8,20,21,24,25]. A simpler and also practical one is 
to choose some specific external parameters such that the 
dynamic phases of the pair of cyclic states accumulated 
in the whole evolution may be eliminated. Interestingly, 
with this method the corresponding cyclic states in NMR 
systems are dark states (the eigenstate with the zero en- 
ergy eigenvalue) , and thus no dynamic phase is involved. 
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The dark state method was proposed for geometric quan- 
tum computation with trapped ions [21], and then de- 
scribed in NMR systems [24] . The other is referred to as 
a two-loop method: let the evolution be dragged by H 
along two closed loops, with one being in t E [0, r] and 
the other in t 6 (r, r + t'\ . The dynamic phases accu- 
mulated in the two loops may be canceled, while the AA 
phases will add. 



III. APPLICATION TO NMR SYSTEMS 

So far, we have proposed a general scheme to achieve a 
universal set of quantum gates based on a pair of orthog- 
onal cyclic states. It is important to further consider im- 
plementing this scheme with real physical systems. Here, 
we illustrate this implementation using NMR systems 
[4,20]. Nevertheless, it is worth pointing out that, in 
principle, the above theory may be applicable to other 
systems which are potentially viable quantum computer 
models. 

For NMR systems, the magnetic field in Eq.(l) or (2) 
in a rotating magnetic field is given by 



B(i) = (B cos cut, Bq amcut, B\), 



(17) 



where £?o,i and u are constants. The qubit state \ip(t)) 
is described by the Schrodinger equation 



(18) 



where the Hamiltonian for a single qubit is given by 

H = -{ljq(t x cosidt + Lu (Ty sincut + uj\o z ) (19) 

with uJi = —g^iBi/h (i = 0, 1) and g being the gyromag- 
netic ratio. The Schrodinger equation with Hamiltonian 
(19) can be solved analytically [23,28]. In terms of ex- 
plicit form of the solution n(x,wi) represented in Ref. 
[23], it is found that a pair of orthogonal initial states 
\ip±) with x — arctan[wo/( w i — w )] take cyclic evolutions 
with the period t = 2it/uj [24], and the evolution paths 
are the curves on a Bloch sphere swept by unit vectors 
±n(x, wt). Therefore, we may use this pair of \if)±) to 
achieve single qubit gates described in Eq.(ll), where 
the corresponding phases for one cycle are given by 



la 

Id 
1 



-7T(1 



UJ 



), 



UJq + wi(wi - lo) 

-7T — . 

LOll 

-7t(1 + Q/w), 



(20) 

(21) 
(22) 



{uio, u>i,u>} satisfying Eq.(12) to accomplish two noncom- 
mutable qubit gates. 

A similar method may be employed to achieve the two- 
qubit operation. The spin-spin interaction in NMR is 
very well approximated by 

Hj = Ja z cr z /2. 

The state of control qubit is (almost) not affected by any 
operation of the target qubit if w\ of the target qubit is 
chosen to be significantly different from uj\ of the control 
qubit. We may prove that the initial states \ip±) de- 
scribed by x S = arctan[cj /(o>f — u))] are a pair of orthog- 
onal cyclic states, and may be used to achieve a two-qubit 
gate described by Eq.(13). Here ui\ = oj\ + (25 — 1) J, ui, 
u>o, and u>\ are parameters for the target qubit (the su- 
perscript "t" is omitted for brevity). The corresponding 
phases for one cycle are given by 



1 S = 
ig 



= -7T(1 



7^ = 
I 5 = 



■«> S M 



-n(l + n s /Lu), 



(23) 

(24) 
(25) 



with n 5 = ^JujI + (w* -uj) 2 . It is seen from Eq.(25) that 
the gate described by Eq.(16) may be accomplished by 
choosing the following special parameters, 



UJ 



UJ = wi 



J. 



It is worth pointing out that we may achieve the nona- 
diabatic geometric gates by choosing some specific pa- 
rameters, with which jd = in the whole process. It is 
direct to verify that the dynamic phase in Eq.(21) is zero 
under the following condition, 



(26) 



Thus the single qubit gates with the parameters satisfy- 
ing Eq.(26) are geometric quantum gates with geometric 
phase shift j g — — 7r(l + w /-\/ w o + ^ ne geometric 
phases 7 S versus the ratio Wi/wo are plotted in Fig. 1. 
It is seen that the nontrivial phases required for two uni- 
versal single-qubit gates may be simply realized by any 
two processes with different values of uji/ljq (except for 
zero or infinite). 

Besides, the geometric two-qubit gates are realized 
whenever 



LU = 2u>i, 

2 2 , t2 



(27) 
(28) 



with f2 = a/wq + (bj\ — uj) 2 . In the derivation of the dy- 
namic phase, E + (t) = \io\ cosx + sinx]?i/2 is used. 
We may choose any two processes with different values 



Correspondingly, the conditional geometric phases are 
given by 7^ = — 7r(l + \Juj{/2uj). Figure 2 shows the 
conditional phases 7^ versus the ratio io\j J. It is evident 
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that the nontrivial phases (7° 7^ 7* in mod 2ir) required 
for two-qubit gates may be achieved for < ui\/ J < 00. 

It is worth pointing out that the constraint described 
by Eq.(26) [ or Eqs.(27) and (28)] is equivalent to the 
condition that the instantaneous dynamic phase for the 
wave function in single qubit (or the target qubit) is al- 
ways zero [24], namely, the states \ip±) used here are the 
dark states. 

The advantage of the above nonadiabatic gates is that 
u> is of the same order of magnitude as loq or u>\. This 
implies that the speed of geometric quantum gate is com- 
parable with that of the dynamic quantum gate. In 
contrast, the speed of quantum gate based on adiabatic 
Berry's phase is much lower than that of gate using dy- 
namic phase, since the adiabatic condition requires that 
both luq and loi should be much larger than u>. 

Note that, the geometric gates U{ \ may not be prac- 
tical by directly using the field B in Eq.(17) as the cor- 
responding geometric phase in Eq.(20) is determined by 
the angle x- For example, Eq.(20) can be rewritten as 
-f g = — 7r(l — cosx); thus -f g = (— 7r) as \ — ( 7r /2)- 
This problem can be solved by rotating the field. It will 
be seen below that the parameter x f° r the initial cyclic 
state may be changed by rotating the symmetric axis of 
field (17), while the phases in Eqs.(20),(21), and (22) are 
invariant. 

We introduce a rotation operator R(y,x' — x) that 
represents the rotation of angle x' — X around the y-axis, 
that is, 



R (2) (y, x' - x) = ex P H(x' - xK] 

in the SU(2) representation, and 

R (3) (y,x' - x) = ex P H(x' - xH 

with 



t 2 




(29) 



(30) 



(31) 



in the 5*0(3) representation. Assuming the required an- 
gle is x' m Eq. (11), we may apply a magnetic field 
B' = R( 3 '(y, x' — x)B, then the solution of the cyclic 
states are n' ± = R( 3 )(y,x' - x)[± n ±(x, 0^±) = 
R-^H^jX' — x)IV'±)] because of the spherical symmetry 
of the system. Thus x ma y change to any required x' for 
implementation of the quantum gate, with the geometric 
phase being unchanged, because the area swept by n' is 
the same as that by n. On the other hand, we have 

£4(t) = -/a4(t)-B'(t) 

= -MR (3) (£,X' - X)n±(t)] • [R (3) (y,x' - x)B(t)} 
= E ± (t). 

Therefore, we have proven that the invariant of all phases 
with respect to the rotation of the symmetric axis of the 



field in Eq. (17). We conclude that 7 (x) in gate (11) is 
determined by the values of {loq,loi,lo} (the symmetric 
axis of the magnetic field). For example, if the magnetic 
field is B' for x' = ( 7r /2) ; we may achieve the geomet- 
ric gate U\ (U2) with 7^2 = — n(l — cosx). It is worth 
pointing out that the above method to control x and 7 g 
separately in quantum gates is also feasible in nongeo- 
metric gates. 

We turn to address how to remove the dynamic phases 
in a multiloop nonadiabatic evolution. The possible gen- 
eralization of a multiloop method from the adiabatic evo- 
lution [8,20] to nonadiabatic case was mentioned in Rcfs. 
[24,25]. We here wish to demonstrate explicitly one re- 
moval procedure of the dynamic phase. 

Let us first choose the magnetic fields in two loops as 

Loop 1 : B= (B coscjt,B smujt,Bi), t £ [0, r) (32) 
Loop 2 : B' = R (3) (j/,c/ - a)(-B' cos Lot, -B' Q smLot, -B[ 



t G [r, 2t] 



(33) 



where r = 2-k/uj, a = arctan[wo/(^i — w)], and af = 
arctan[wo/(w'! + co)] with w\ = —g\iB' i j'h (i = 0,1). As 
shown before, a pair of orthogonal initial states \ip±) 
(IV'i)) w ith x = a (x' = a ') take cyclic evolutions dur- 
ing the loop one (two). The rotation R^(y, a' — a) in 
Eq.(33) ensures that the cyclic initial states in the two 
loops are the same at the time t — 2it/uj [31]. There- 
fore, the gate described by the two loops is given by 
U = U{x,l {1) + 1 {2) ), where 7 (1) (7 (2) ) is the total 
phase accumulated in the loop one (two). Denoting 7^ 

{I — 1,2) and 7^ the dynamic phases and geometric 
phases accumulated in the loop I, respectively, we now 
illustrate that there exist processes satisfying 



-IV 



0. 



where — Tir is a nontrivial geometric phase which we in- 
tend to realize in geometric quantum gates. Then the 
magnetic fields should satisfy the following equations 



LOl — LO U>[ + LU 



where 



ft 



ft' 



2-T, 



K) 2 + K) 2 



■w 



(34) 
(35) 
and 



ft' = y/ (i^o) 2 + ( w 'i + UJ ) 2 - As for the required T, it is 
possible that there exist many solutions, since there are 
five unknown variables in two equations. For example, 
we numerically calculate the solutions for T = 1/2. For 
simplicity, we set oj[ = uii as the unit, and find that if 
{w, uj ,uj' } (0 < uj < LOi) satisfy the equations given by 

to + 1.13389wo = 0.99998, 

Lo + 1.07091w - 0.06299^0 = 0.88889, 



G 



which describe a straight line(segment) in the three- 
dimensional space, the geometric phase accumulated in 
the whole two-loop evolution is just what we required, 
with the total dynamic phase being zero. 

The multiloop method to remove dynamic phase is 
also feasible for two-qubit geometric quantum gates. We 
choose the magnetic fields on the target qubit in two 
loops as 

Loop 1 : B= (B coswt, -Bo sinwi, B\), te[0,r), (36) 
Loop 2 : B' = R^\y,7])(-B' cosw't, -B' Q smu>'t, -B[), 
te[r,r + 4 (37) 

where r' = 2w/u)'. The angle 

r\ = arctan[w /(w 1 — u>)] — arctan[w / (co'-f + uu')] (38) 

should be independent on the state 5 in the control qubit. 
To guarantee that the interaction between qubits is still 
determined by the original initial state 6 of the control 
qubit, the control qubit should be rotated by R^ 3 \y,r]) 
at time t = r [The state of the controlled qubit is un- 
changed if a rotation R^(y, —77) is also applied at time 
t + t']. Correspondingly, the S— independent constraint 
described by Eq.(38) can be rewritten as 



L0 



{iOr-uof-J 2 (u)[ + uu') 2 - J 2 ' 



(39) 



On the other hand, the condition under which there exist 
processes with zero dynamic phase is 



lu 2 + u{ {uj{ - u) _ K) 2 + Lj[ d + J) 



to 



(40) 



where fl s = V^q + ( w i — w ) 2 an d 

fl' 5 = \/(uj' ) 2 + {uj[ s + u') 2 . Note that the geometric 
phases are nontrivial ( 7<J 7^ 7° in mod 2tt), and thus can 
be applicable in geometric quantum computation. 

The magnetic fields, which satisfy Eqs. (39) and (40) 
in loop two, as a function of u are plotted in Fig. 3a, 
where lo = loi = 5.0 with J as the unit. We may numer- 
ically calculate the three unknown variables {w',cj ,wi} 
in three equations described by Eqs. (39) and (40). Then 
the conditional geometric phases may be obtained from 
equations 7* = — T 5 n with 



The corresponding conditional phases 7* for uj = w\ — 
5.0 as a function of to are plotted in Fig.3(b). It is 
seen that the nontrivial phases 7^ ^ 7° may be realized 
by appropriately choosing the values of {iv 7 ujq,uji} and 
{iv 1 , iv'q, As a consequence the nontrivial two-qubit 
geometric quantum gate may be achieved. 



IV. CONCLUSIONS AND DISCUSSIONS 

An experimentally feasible scheme based on a pair of 
orthogonal cyclic states has been proposed to accomplish 
a universal set of quantum logical gates, in which quan- 
tum computation implemented by both dynamic and ge- 
ometric operations can work, i.e., quantum gates in this 
scheme may be implemented by the total phases accumu- 
lated in the cyclic evolution, and the geometric quantum 
computation can be achieved by eliminating the dynamic 
phase. Furthermore, the geometric phase shift used is 
the cyclic AA phase, which can be nonadiabatic. It is 
possible that the gates achieved here can handle arbi- 
trary quantum computation without the intrinsic limi- 
tation on operation time. Therefore, the nonadiabatic 
method proposed here may allow us to physically imple- 
ment (geometric) quantum computation even for systems 
with very short decoherence time, which could be espe- 
cially useful for solid-state implementations of scalable 
quantum computers. 

We here discuss briefly the errors induced by ran- 
dom noises in geometric quantum computation. Random 
noises may lead to two kinds of errors. One is that the 
path may not be exactly closed at the end of the gate op- 
eration, leading to the noncyclic corrections. The other 
is that the evolution path may fluctuate around the ideal 
path with known cone angle. The noncyclic corrections 
could be negligible at least when the first-order correc- 
tions from random noises are taken into account, as was 
indicated in Ref. [32]. On the other hand, as that in the 
adiabatic cyclic geometric scheme, the present scheme 
is also robust to the second type of errors as the area 
enclosed by the evolution path (geometric phase) is in- 
sensitive to the random fluctuation. 

Finally, we wish to make a few remarks on experimen- 
tal implementation of geometric quantum computation. 
The simplest geometric quantum computation should ex- 
perimentally complete the following three steps one by 
one: (i) detection of the (conditional) geometric phase 
shifts in qubit systems; (ii) implementation of a univer- 
sal set of geometric quantum logic gates, particularly the 
implementation of a conditional gate, (iii) illustration 
of a simple algorithm by pure geometric quantum gates, 
such as Deutsch's problem, Grover's search algorithm, or 
Shor's factorization algorithm, etc.. Two recent excit- 
ing experiments reported that the conditional geometric 
phase shifts for quantum logical gates using NMR were 
detected in adiabatic [20] and nonadiabatic [2G] regions. 
However, a universal set of gates as well as a simple quan- 
tum algorithm experimented by (adiabatic or nonadia- 
batic) geometric phases are still awaited. 
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FIG. 1. The geometric phase 7 3 versus the ratio wi/wq. 
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FIG. 2. The conditional phases 7* versus the ratio u>i/J. 



FIG. 3. (a) The magnetic fields required in loop two versus 
ui. (b) The conditional geometric phases 7^ versus uj. 
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